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^ ' ABSTRACT 

D ■ Context. The variable Sun is the most likely candidate for natural forcing of past climate change on time scales of 50 to 1000 years. 

Evidence for this understanding is that the terrestrial climate correlates positively with solar activity. During the past 10000 years, the 
Sun has experienced substantial variations in activity and there have been numerous attempts to reconstruct solar irradiance. While 
there is general agreement on how solar forcing varied during the last several hundred years — all reconstructions are proportional to 
the solar activity — there is scientific controversy on the magnitude of solar forcing. 

Aims. We present a reconstruction of the Total and Spectral Solar Irradiance covering 130 nm-10 yum from 1610 to the present with 
\ annual resolution and for the Holocene with 22-year resolution. 

Methods. We assume that the minimum state of the quiet Sun in time corresponds to the observed quietest area on the present Sun. 
Then we use available long-term proxies of the solar activity, which are 10 Be isotope concentrations in ice cores and 22-year smoothed 
neutron monitor data, to interpolate between the present quiet Sun and the minimum state of the quiet Sun. This determines the long- 
O { term trend in the solar variability which is then superposed with the 1 1 -year activity cycle calculated from the sunspot number. The 

time-dependent solar spectral irradiance from about 7000 BC to the present is then derived using a state-of-the-art radiation code. 
Results. We derive a total and spectral solar irradiance that was substantially lower during the Maunder minimum than observed 
today. The difference is remarkably larger than other estimations published in the recent literature. The magnitude of the solar UV 
variability, which indirectly affects climate is also found to exceed previous estimates. We discuss in details the assumptions which 
leaded us to this conclusion. 

Key words. Sun: solar-terrestrial relations - Sun: UV radiation - Sun: atmosphere - Radiative transfer - Line: formation -Sun: 
surface magnetism 
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' 1 . Introduction question of construction of the unique self-consistent TSI corn- 

el" . _ _ . , , , . posite for the satellite epoch is still open. 

. The Sun is a variable star whose activity vanes over time- »„ ^ , n i , T j- 

i c ■ . . ii • r\ .u i . .u- . Measurements of the Spectral Solar Irradiance (SSI) are even 

Cn scales ranging from minutes to millennia. Over the last thirty , , , 

O., , -j. , , more difficult and instrumental problems prevent the construc- 
years the solar irradiance was measured by numerous space - . . j T , . r * , . ,im = r~n 

— — I > . . . c rr . , c , T j. /rrcifu tion of essential composites Kxivova et al. 2011; Domingo et al. 

missions. Measurements of the Total Solar Irradiance (TSI) be- u»»» l „ , rrr-P , k i — — — : — , „ ™ mr 

t-H ., . , , u f .u xmmnc 7 ■ ■ • 2009). Recently dKnvova et al. 2009b) the theoretical SATIRE 

. . ■ came available with the launch of the NIMBUS 7 mission in ~z — , . , J " , T trr — , arr~- n 

. lr no Jtt — t > i ll ir.r.o K o- .u t, ct (Spectral And Total Irradiance REconstruction) (Krivova et al. 

> 1978 Hoytet al. 1992). Since then TSI was measured by sev- \ n ^ rr-. „ c , — r. OAAOi j, . CITCI u /c i 

, : . . . „ , c ., a At a 2003; Krivova & Solanki 2008) model and SUSIM (Solar 

eral consecutive instruments. Each of them suffered from degra- 77; — . , — — ,. . . . 

, . , . ,. . , , rr • j . Ultraviolet Spectral Irradiance Monitor) measurements were 

dation and individual systematic effects, so the direct compar- , r , , . . . . . ,. , , 

$— i . „ , ... ... rp, TCT .. used to reconstruct the solar UV irradiance back to 1974. 



, , 2003; Krivova & Solanki 2008) model and SUSIM (Solar 

* j aril pAnrociitnfo incrnimonrc h o/<h (it thnm en rro ra ri Tr/\ni rtorrm- v 

^3 ' ison of the measurements is impossible. Three TSI composite lu iwwmuuM ^ w ' •' '"■^•"■"^y ' v / _r ' 

based on the available data were constructed by three groups: Takln 8 lnt0 account the Problems of the reconstruction of 

PMOD (Frohlich 2006), ACRIM (Willson & Mordvinov 2003), TSI and SSI dunn 8 the recent P enod of satelllte observations 

and IRMB (Dewit te^ECKMh . These composites give quite and combining it with the fact that the solar dynamo, which is 

different values of TSI, especially before 1980 and during the belleved to drive all activity manifestations, is still not fully un- 

so-called ACRIM gap between June 1989 (end of ACRIM I derstood ( |Charbonneau| | 2010|) , one recognizes the difficulty of 

observations) and October 1991 (beginning of ACRIM II ob- ^constructing TSI and SSI to the past, when no direct measure- 

servations). The most striking detail is the increase of TSI be- ments were avallable - 

tween the minima 1986 and 1996 in the ACRIM composite and Long-term changes in solar irradiance were suspected as 

the absence of such increase in the PMOD and IRMB compos- earl y as the mid-nineteenth century (|C. P. Smyfh|[1855|). One 

ites. Although significant progress was made during the last few of the first quantitative estimates of its magnitude as well as 

years and the increa se of the TSI in the ACRIM composi te was P ast solar irradiance reconstructions was obtained by using the 

strongly criticized (Frohlich 2009; K rivova et all 12009a) , the observations of solar-like stars (|Lean et al.||l995|). It was con- 

eluded that TSI during the Maunder minimum was about 3- 
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diative forcing] A,Fp_M~ 0.5-0.7 W/m 2 . However, these results 
were not confirmed by large survey s of solar-like stars and ar e 
no longer considered to be correct lHall & Lockwoodl d2004). 
Recent re constructions based on the magnetic fi eld surface dis- 
tribution (Wan g et al.ll2005l; iKrivova et alj|2007l) and on extrap- 
olation of the assumed correlation b etween TSI and the open 
magn etic flux (|Lockwood et alJll999h durin g the last three min- 
ima (Steinhilber et al. 2009; Frohlich 2009) resulted in a low so- 
lar forcing value within the range AF P _ M ~ 0.1-0.2 W/m 2 . 



2. Effects of solar radiative forcing on climate 

Variations on time-scales up to the 27-day rotational period have 
an important influence on space weather but not on terrestrial cli- 
mate. The effects of the 1 1 -year solar cycle are clearly detected 
in the atmosphere and are widely discussed in the literature (e.g., 
Egorova et al. 2004; Haigh 2007), while the imp rints from vari - 
ations on longer time-scales are more subtle dGrav et al.ll2010h . 
Analysis of historical data suggests a strong correlation between 
solar activity and natural climate variations on centennial time- 
scales, such as the colder climate during the Maunder (about 
1650-1700 AD) and Dalton (about 1800-1820 AD) minima as 
well as climate warming during the steady increa se in solar ac- 
tivity in the first half of the twentieth-century (|Siscoel 1978t 
iHovt & Schatte73fl997t ISolomon et al]l2007t iGrav et aUboiOl) . 
Numerous attempts to confirm these correlations based on dif- 
ferent climate models have shown that it is only possible if ei- 
ther the applied perturbations of direct solar radiative forcing 
are large (consistent with a direct solar radiative forcing from 
the present to Maunder minimum AFp_M~ 0.6-0.8 W/m 2 ) or 
the amplification of a weak direct solar forcing is substantial. 
Because the majority of recent Afp_M estimates (see Sect. 1) 
are only in the range 0.1-0.2 W/m 2 , and amplification processes 
have not been identified, the role of so lar forcing in natura l cli- 
mate change remains highly uncertain (Solomon et al.l2007l) . In 
this paper we show that the solar forcing may be significantly 
larger than reported in the recent publications. 



3. Methods 

In this paper we present a new alternative technique, avoiding 
calibration of our model with presently observed TSI variations 
and extrapolation to the past. We assert that the amou nt of mag- 
netic energy that remains present (Ide Wiin et al.l20 09) at the sur- 
face of a spotless (i.e. quiet) Sun is the main driver of solar ir- 
radiance variability on centennial time scales. The main concept 
of our technique is to determine the level of the magnetically en- 
hanced contribution to the irradiance of the present quiet Sun. 
Then, for the reconstruction to the past, this magnetically en- 
hanced component has to be scaled with the proxies for the quiet 
Sun activity. However a proxy for the long-term activity of the 
quiet Sun does not yet exist. A good candidate for such a proxy is 
the sm all-scale turbu lent magnetic fields accessed with the Hanle 
effect dStenfloll9 82). However the consecuti ve measurements o f 
these fields are limited to the last few years (iKleint et al.ll2010b . 
Therefore, we assume that the existing proxies of the solar activ- 
ity averaged over the two solar cycles period can also describe 
the activity of the quiet Sun. Averaging of proxies allows suf- 
ficient time for the magnetic components to decay to the quiet 



1 Solar radiative forcing is a direct energy source to the Earth and is 
related to the change in TSI by AF = ATSI ■ (1 - A)/4, where A is the 
Earth's albedo. 



network, and then to even smaller mag netic features as the de- 
cay process can take up to several years (ISolanki et al.ll2000l) . In 
other words we set the small-scale activity (which defines the 
fractional contributions of quiet Sun components) to be propor- 
tional to the large-scale activity. We want to state clearly that this 
proposition is an assumption, which is however in line with a 
high-resolution observations of the present Sun. The large-scale 
structure due to strong magnetically active features is repeated 
on a smaller scale in less active regions, and even in the ap- 
parently quietest areas there is still a mosaic of regions of dif- 
ferent magnetic field strengths, reminiscent of fractal structure 
(Ide Wiin et al. l l2009h . 

Hence, the time-dependent irradiance / qu iet(^, of the quiet 
Sun in our reconstruction can be calculated as 



In 



1 quiet 



(A, to) - / m in.state(T) 



<Proxy > 22 (t) 



<Proxy > 22 (to) 



(1) 



where / qu i e t(/i, t) is the time dependent irradiance of the quiet Sun 
and to denotes a reference time. / Proxy 11 is the averaged over 
a 22-year period value of the proxy for the solar activity. The 
22-year period was chosen because the cosmogenic isotope data 
used for the re construction (see below ) are available as 22-year 
averaged data ( Steinhilber et al. 2008). We have set 1996 as the 
reference year, i.e. 7 qu ietW, to) is the irradiance of the quiet Sun 
as observed during the 1996 minimum. Let us notice that the 
quiet Sun irradiance was roughly constant for the last 3 cycles 
(see below). Therefore the solar spectrum during the 1996 mini- 
mum is a good representation of the present quiet Sun spectrum. 
/n J in.state( / i s me irradiance of an absolute minimum state of the 
Sun, with a minimum of remaining magnetic flux emerging on 
the solar surface. Thus, the term in brackets can be considered 
as the enhancement level of the present quiet Sun with respect to 
the Sun in its most inactive state. 

The prominent, readily observable active regions on the Sun 
also contribute to the variability in irradiance. Thus, the full solar 
variability is described by: 



I(A, t) = / qu ie t (i, t) + /activeW, t), 



(2) 



where I ac &ve(A, t) is the contribution to solar irradiance from ac- 
tive regions e.g. sunspots, plages, and network. /active^ , i s cal ~ 
culated following the approach bv lKrivova et alJ (l2003h (see also 
the online Sect. [6} and, as calculations are done with annu- 
lar resolution, is proportional to the sunspot number. The group 
sunspot number used in our r econstruction is taken f rom NOAA 
data center and described in IHovt & Schattenl (119981) . 

The term /active^, in the Eq. (2) describes the cyclic com- 
ponent of solar variability due to the 1 1 -year activity cycle, while 
the slower long-term changes due to the evolution of the small- 
scale magnetic flux are given by the Eq. (1). 4ctiveW, term can 
only be calculated for the time when sunspot number is avail- 
able and therefore our reconstruction has 22-year resolution for 
the Holocene and annual resolution from 1610 to present. 

The choice of the model for the minimum state of the Sun 
is a crucial point in our technique as it defines the amplitude 
of the reconstructed solar irradiance variability. Observations 
of the Sun with relatively high spatial resolution show that the 
present quiet Sun is still highly inhomogeneous even though 
it appears spotless. Measurements obtained with the Harvard 
spectroheliometer aboard Skylab were used to derive bright- 
ness components of the quiet Sun and to constr uct correspond- 
ing semi-empirical solar at mosphere structures (Vernazz a et alJ 
1 198 It iFontenla et al.|[l999l) . The darkest regions with the least 
amount of magnetic flux corresponds to the faint supergranule 
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cell interior (component A). This component comes closest to 
describing the most inactive state of the Sun. Therefore, in our 
approach we set / m i n . s tateW) = Ia(X)- Thus, the basic magnitude 
of solar irradiance variations is given by the difference between 
the irradiance of the present quiet Sun (composed from a dis- 
tribution of brightness components defined in supporting online 
material) and the irradiance from component A (see Eq. (1)). 

For the reconstruction to the past this amplitude is scaled 
with proxies for solar activity. Two proxies are available for 
the reconstruction: Group sunspot number, which is available 
from the present to 1610 AD, and the solar modulation poten- 
tial extending back to circa 7300 BC. The latter is a measure 
of the heliospheric shielding from cosmic rays derived from the 
analysis of cosmogenic isotope abundances in tree rings or ice 
cores , and is available with a time resolution of 2-3 solar cy- 
cles dSteinhilber et al.ll2.008h . Although sunspot number dropped 
to zero for a long time d uring the Maunder minimum, the solar 
cycle was uninterrupted (Be er et al.l ll 998b lUsoskin et al.l l2001 ) 
and the modulation potential did not fall to zero. Hence, a recon- 
struction based solely on sunspot number may underestimate the 
solar activity during the Maunder minimum. Therefore in our re- 
construction we used the solar modulation potential to calculate 
the long-term variations and sunspot number to superpose them 
with the 1 1 -year cycle variations (see the Online Section 6.2). 

The modulation potential used in the calculations is based 
on the composite of data determined from the cosmogenic iso- 
tope records of 1() Be and neutron monitor. 10 Be data are available 
up to about 1970 (McCrack en et al.l 12004) and neutron moni- 
tor data, which are used to calculate the current solar modu- 
lation potential, are available since the 1950s. Three different 
datasets of the 10 Be data were used: measure ments at DYE 3, 
Green land, and the South Pole provided by iMcCracken et al.1 
(2004j) for the reconstruction back to the Maunder minimum, 
an d measurements from th e Greenland Ice core Project provided 
by IVonmoos et al.l d2006h for the reconstruction back to 7300 
BC. T he neutron monitor data were provided by Usoskin et al. 
(120051) . Although they are available with a monthly resolution, 
we calculated the 22-year mean of the neutron monitor in order 
to homogenize the data sets and because only averaged modula- 
tion potential can be used as a proxy for the quiet Sun activity. 
The transition from 10 Be to neutron monitor data is shown in the 
lower panel of Fig. 1 . Prior to 1952 the composites are based on 
10 Be data (cyan based on DYE 3, and red on South Pole records, 
both with 22-year resolution). Values between 1952 and 1998 are 
obtained by linearly interpolating the following consecutive data 
points: two blue crosses at 1952 (the last data points used from 
10 Be records), green crosses at 1976 and 1998 (22-year averages 
of the neutron monitor data). Values after 1998 require knowl- 
edge of the next 22-year average (between 2010 and 2031) of 
the modulation potential. We assume this average to be 92% of 
the previous average (between 1988 and 2009), which allows us 
to reproduce the observed TSI minimum in 2008. The datasets 
mentioned above were derived using different assumptions of 
the Local Interstellar Spectra (LIS). To homo genize the data 



we co nvert ed the modula ti on po tential from McCracken et al. 



(20041) and lUsoskin et ail (120051) to LIS bv ICastagnoli & Lai 
(1980S) which is used in IVonmoos et al.l d2006l). For conver - 
sion we appl i ed a m ethod suggested bv lSteinhilber et a"D d2008). 
iHerbst et al. (2010) analyzed the dependency of the modulation 
potential on the applied LIS models and show that negativ e val- 
ues of the modulation potential in Vonmo os et alJ J2006) data 
can be corrected with another LIS model. Let us notice that the 
solar forcing in our reconstruction is determined by the rela- 
tive values of the modulation potential so change to the differ- 
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Fig. 1. Modulation potential and TSI reconstruction for the last 
70 years. Lower panel: Yearly averaged neutron monitor data 
(green) and two modulation potential composites (red and cyan 
curves, see the discussion in the text). Upper panel: TSI recon- 
structions based on the two modulation potential composites (red 
and cyan curves). The black dashed line is the observed TSI from 
the PMOD composite. Error bar for 1980 corresponds to 1/4 the 
difference between two published TSI c omposites, and the error 
for 2008 is taken from iFrohlichl (12009). The reconstructed TSI 
curves are normalized to the 1996 minimum and the grey-shaded 
region indicates the intrinsic uncertainty due to differences in the 
modulation potential data. 



ent LIS model will introduce only several percents correction, 
which is much less than estimated accuracy of the reconstruc- 
tion. In our dataset the value of the reference modulation po- 
tential <Proxy >2i(to) (see Eq. (1)) is equal to 631 MeV The 
relativ e error is estimated to be less than 10% (lUsoskin et al .1 
120051) . 

Synthetic solar spectr a are calculated with a state-of-the- 
art ra diative transfer code (lHaberreiter et alJ 2008 ; Shapir o et al.l 
2010). We discuss the calculations in more detail in online 
Sect. [6] where in addition we show that the concept of quiet 
Sun activity scaling can also be expressed in terms of varying 
fractional contributions from different components of the quiet 
Sun. 



4. Results and Discussion 

In the upper panels of Fig. 1 and Fig. 2 we present the TSI re- 
constructions, which are obtained after the integration of the 
Eqs. (1) and (2) over the wavelengths and normalization of the 
quiet Sun value for the reference year 1996 to 1365.5 W/m 2 . As 
the sunspot number is only available since 1610 AD the recon- 
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Fig. 2. Modulation potential (lower panel) and TSI reconstructions (upper panel) for the last 2500 years. Data prior to 1600 AD are 
based on the modulation potential derived from 10 Be records from the Greenland Ice core Project (red curves). Data since 1600 AD 
are based on the two composites shown in Fig. 1 (red and cyan curves). The grey-shaded area indicates the intrinsic uncertainty. 



struction of the full solar cycle variability with an annual reso- 
lution extends back only 400 years. Both reconstructions in the 
right-hand panels of Fig. 2 are based on the 10 Be data sets men- 
tioned above. The difference in the reconstructions allows the 
error originating from the uncertainties in the proxy data to be 
estimated (20-50 % in the solar forcing value, depending on 
the year). This is large but still significantly less than the change 
in irradiance between the present and the Maunder minimum. 
Both reconstructions suggest a significant increase in TSI dur- 
ing the first half of the twentieth-century as well as low solar ir- 
radiance during the Maunder and Dalton minima. The difference 
between the current and reconstructed TSI during the Maunder 
minimum is about 6 + 3 W/m 2 (equivalent to a solar forcing of 
AFp_M~ 1.0+0.5 W/m 2 ) which is substantially larger than recent 
estimates (see Sect. [TJ. Note that as our technique uses 22-year 
means of the solar modulation potential our approach cannot be 
tested with the last, unusual solar minimum in 2008. In order 
to reproduce the current minimum as shown in Fig. 1 we have 
adopted a value of 584 MeV for the future 22-year average in 
2020 (which is 92% of the 22-year average for 1988-2009). 

The reconstruction before 1500 AD in the left-hand panel 
of Fig. 2 (which stops for clarity at 500 BC) is based on 10 Be 
record s from the Greenland Ice core Project dVonmoos et al.l 
2006). The modulation potential during the Maunder minimum 
is about 3^1 times less than at present but not zero. However, 
it has decreased to zero several times in the past and the corre- 
sponding TSI was even smaller than during the Maunder min- 
imum. There were also several periods when the modulation 
potential, and hence TSI, were higher than the present value. 
The reconstruction back to 7000 BC is presented in Fig. 3. The 
choice of model A introduces an uncertainty of the order of 30%, 
which is estimated by comparing model A to other possible can- 
didat es for the min imum state of the quiet Sun, e.g., model B 
from Vernazza et al.1 (Il98ll) . Combining this with the uncertain- 
ties of the proxy data outlined in Fig. 2 we can roughly esti- 
mate the uncertainty of our solar forcing value to be 50%. In 
additional to the 10 Be-based data of solar activity there are sev- 
eral 14 C-based datasets (e.g., ISolanki et alJl2Qp4tlVonmoos et all 
2006; iMuscheler etafl |2007t |Usoskin| |2008|). Employment of 
these datasets will lead to a somewhat different values of the 



solar variability, which is, however, covered by our rough order 
of magnitude estimate of the overall uncertainty of the recon- 
struction. 

Our TSI reconstructions give a value of ~1 W/m 2 per decade 
for the period 1900-1950. The Smithsonian Astrophysical 
Observatory (SAO) has a 32-year record of ground-based ob- 
servations fo r 1920-195 2. Although the SAO data are disputed 
in reliability (Abbo tl2~007l) and clearly contain a non-solar signal 
they are the only available long-term measurements of TSI in 
the first half of the twentieth-century. The data show an increase 
of 1+0.5 (1.5+0.5) W/m 2 per decade for the period 1928-1947 
(1920-1952). As we are aware that maintaining a stable calibra- 
tion to better than 0.1 % over 30 years is very demanding we 
cannot claim that the historical data confirm our reconstruction. 
Nevertheless, it is intriguing to note how well the SAO trend 
agrees with our TSI reconstruction. 

Our reconstructed solar spectral irradi- 
ance comprises spectra from 130 nm to 
10 /jm. Fig. 4 presents a reconstruction of the integrated 
flux for several, selected spectral regions. The contrast between 
different brightness components of the quiet Sun is especially 
high in the UV, which results in a large historical variability of 
the UV spectral irradiance. The irradiance in the Schumann- 
Runge bands and Herzberg continuum increases from the 
Maunder minimum to the present by about 26.6% and 10.9% 
respectively, which is much larger than 0.4% for TSI and the 
visible region. The variability is also relatively high around the 
CN violet system whose strength is very sensitive to even small 
temperature differences due to the high value of the dissociation 
potential. The large UV variability reported here is especially 
of importance to the climate community because it influ ences 
climate via an indirect, n on-linearly amplified forcing (Haigh 
Il994t lEgorova et al.l[2004h . 

We are aware that the choice of model A is responsible 
for a relatively large fraction of the uncertainty in our results. 
Higher resolution observations have recently become available, 
and model A could possibly be improved in future studies. We 
emphasize that model A is not the coldest possible quiet Sun 
model and therefore our estimate is not a lower limit of the Sun's 
energy output. The coldest model would be a non-magnetic at- 
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Fig. 4. Reconstruction of the integrated spectral irradiance in se- 
lected wavelength bands. This reconstruction is based on the 
composite of the I0 Be South Pole record and neutron monitor 
data (red curve in Fig. 2). Panels from bottom to top: The 500- 
600 nm band is representative of the time evolution of the visi- 
ble irradiance. The CN violet system is one of the most variable 
bands accessible from the ground. The Herzberg continuum is 
crucial for ozone production in the stratosphere. The Schumann- 
Runge band is important for heating processes in the middle at- 
mosphere. 



mospheric structure without a chromosphere and corona. As 
model A contains some remaining magnetic activity our ap- 
proach does not imply that the solar dynamo stops during the 
periods when the modulation potential is equal to zero. Let us 
also notice that the modulation potential never reached zero for 
the last 400 years (see Fig. 2). 



5. Conclusions 

We present a new technique to reconstruct total and spectral so- 
lar irradiance over the Holocene. We obtained a large historical 
solar forcing between the Maunder minimum and the present, as 
well as a significant increase in solar irradiance in the first half of 
the twentieth-century. Our value of the historical solar forcing is 



remarkably larger than other estimations published in the recent 
literature. 

We note that our conclusions can not be tested on the basis 
of the last 30 years of solar observations because, according to 
the proxy data, the Sun was in a maximum plato state in its long- 
term evolution. All recently published reconstructions agree well 
during the satellite observational period and diverge only in the 
past. This implies that observational data do not allow to select 
and favor one of the proposed reconstructions. Therefore, until 
new evidence become available we are in a situation that differ- 
ent approaches and hypothesis yield different solar forcing val- 
ues. Our result allows the climate community to evaluate the full 
range of present uncertainty in solar forcing. 

The full dataset of the solar spectral irradiance back to 7000 
BC is available upon request. 
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Fig. 3. TSI reconstruction from 7000 BC to 500 AD. The reconstruction is based on the modulation potential derived from 10 Be 
records from the Greenland Ice core Project. Red dotted lines indicate the estimated error bars of the reconstruction. Black dashed 
and dotted lines indicate TSI for the 1996 solar minimum and the lowest TSI during the Maunder minimum, respectively. 
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6. Spectral synthesis 

6.1. Quiet Sun 



The quiet Sun is a combination of different brightness compo- 
nents and the evolution of their fractional contributions drives its 
activity and long-term irrad iance variability. The four main com- 
ponen ts of the quiet Sun (IVernazza et al.l [19811 iFontenla et al.l 
1999) are: component A (faint supergranule cell interior), com- 
ponent C (average supergranule cell interior), component E (av- 
erage network or quiet network), and component F (bright net- 
work). 

We calculated the synthetic spectra I a, Ic, Ie, Ip of all these 
components employing the NLTE (non-local t hermodynamic 
equilibrium) COde f or Solar Irradiance (Haberreiter et al. 2008; 
IShapiro et all 1201 Ol) (COSI). Recently, Shapiro et al. (2010) 
showed that COSI calculations with the atmosphere model for 
component C reproduces spectral irradiance measurements from 
the last two solar minima with good accuracy. This is used in Eq. 
(1) from the main text, where /p reS ent(^) is substituted by the ir- 
radiance Ic(A) for component C. 

The evolution of the magnetic activity of the quiet Sun can 
be represented by the time-dependent fractional contributions 
(i.e. filling factors) of different components of the quiet Sun. The 
quiet Sun can be described by model A with varying contribu- 
tions from the brighter components. We compose the quiet Sun 
using model A, and model E which is an adequate representation 
of the brighter contributions. 

Using COSI we demonstrate (see Fig. 5) that the solar irra- 
diance for model C and, therefore, measured solar irradiance for 
the last two minima, can be successfully reproduced with a com- 
bination of advent = 43% model A and o^^t = 57% model E 
(a^ + a E = 1 V 

V^present present 

= IC(A) = Resent W) + <e.sen, (3) 

The linear scaling of the magnetic activity of the quiet Sun 
with proxies, as described in the main text, is equivalent to set- 
ting the model E filling factor to be proportional to a chosen 
proxy. So that the time-dependent irradiance of the quiet Sun 
can be calculated: 



/quietW, = (l - a E (t))l A (A) + a E (t)I E (A), 



where 

a E (t) = ? rOXy{t) a E 



Proxy 



present 



present • 



(4) 



(5) 



The set of Eqs. (2-4) is equivalent to Eq. (1). 

There is an ongoing discussion of whether the trend in the 
filling factor of the quiet netwo rk can be detected. It was sug- 
gested (Foukal & Milano 2001) that the analysis of the histori- 
cal Mt. Wilson observations exclude the existence of the trend 
between 1914 and 1996. H owever, this analysis was criticized 
for using uncalibrated data (Solanki & Krivova 2004) and also 
contradicts other studies (lLockwoodl20 1 Oh . Thus only the future 
long-term monitori ng of the quiet Su n with high resolution or 
recently proposed dKleint et al.l I2010D monitoring of weak tur- 
bulent magnetic fields can help to clarify this question. 
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Fig. 5. Representation of the quiet Sun by a combination of mod- 
els. Ratios of synthetic solar spectra of models A (blue) and E 
(red) to model C. The black line results by combining 57% of the 
model E spectrum with 43% of the model A spectrum. The spec- 
tral signatures are almost perfectly cancelled, implying that the 
combination of both spectra is equal to the model C spectrum. 
The latter has been demonstrated to accurately rep roduce the ob- 
serve d solar spectra during the last two minima( Shapiro et al. 
2010). Hence, a combination of spectra for models A and E re- 
produces the observed quiet Sun spectrum. 



we follow the approach by iKrivova et al. (2003). Because our 
main goal is to reproduce centennial solar variability and be- 
cause magnetograms are unavailable for historical time peri- 
ods, we scale the faculae and the active network filling factors 
with the sunspot number instead of using filling factors derived 
from available magnetogram data. The synthetic spectra are then 
added according to their filling factors, and TSI is then deter- 
mined by integrating over all wavelengths. 



6.2. Active Sun 

From 1610 onward we have additional information from sunspot 
number, which allows the calculation of the active regions con- 
tribution to the solar irradiance (7 ac tive(^ in Eq. (2)). For this 



